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A Method  for  Treating  the  Sheath  Size  in  the 
Langmuir  Mott-Smith  Equations 


1.  INTRODUCTION 

Recent  work  in  probe  theory1  has  been  largely  oriented  towards  collection  of 
ions  from  laboratory  plasmas  in  which  the  electron  temperature  greatly  exceeds 
ion  temperatures.  In  the  ionosphere,  electron  temperatures  are  at  most  only  2 
to  3 times  ion  temperatures  and  are  often  even  lower  than  this.  In  these  circum- 
stances  the  original  Langmuir  Mott-Smith  (LMS)  equations  are  applicable  not  only 
to  electrons  but  also  possibly  to  ions,  and  if  some  means  of  treating  a sheath  of 
unknown  size  were  available,  the  more  recent  (and  more  complex)  theories  with 
the  attendant  difficulties  in  use  could  be  dispensed  with.  This  is  the  problem  which 
is  addressed  in  the  present  paper.  The  treatment  developed  here  will  also  be 
applicable  to  the  charging  of  vehicles  by  emission  of  beams  of  charged  particles 
into  the  ionosphere,  since  with  the  interchange  of  roles  of  dependent  and  independent 
variables,  the  current -voltage  relationships  in  this  case  may  be  identical  with  those 
governing  the  behavior  of  electric  probes. 


(Received  for  publication  2 June  1978) 

1.  Swift,  J.D.  and  Schwar,  M.J.R.  (1969)  Electrical  Probes  for  Plasma  Diagnos- 

tics, American  Elsevier,  pp  69-72. 

2.  Suits,  C.G. , Ed.  (1961)  The  Collected  Works  of  Irving  Langmuir,  Vol.  4, 

Per  gam  on  Press. 
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Although  the  LMS  formulas  have  been  discussed  extensively  in  the  literature, 
there  still  seems  to  be  some  question  concerning  both  the  limits  of  their  applica- 
bility and  also  their  relationship  to  the  Langmuir  space  charge  limited  (SCL) 

O 

equations  for  diodes.  We  are  not  referring  to  limitations  such  as  the  shape  of 
potential  curves  or  relative  values  of  positive  and  negative  particle  temperatures, 
but  rather  to  more  obvious  scaling  effects.  In  much  previous  work,  the  limiting 
cases  of  the  LMS  equations  are  discussed  solely  in  terms  of  the  variation  in  sheath 
size.  Here,  the  limiting  cases  are  discussed  in  terms  of  combinations  of  sheath 
size  and  potential,  as  they  actually  appear  in  the  LMS  equations.  As  a result  of 
this  treatment,  it  will  appear  that  the  relationship  between  the  LMS  and  SCL  equa- 
tions is  not  necessarily  one  of  mutual  exclusion  but  rather,  often,  of  complemen- 
tarity as  far  as  domains  of  applicability  are  concerned.  From  the  two  equations, 
the  sheath  size  can  then  be  eliminated  and  an  expression  obtained  for  the  current 
collected  in  terms  of  voltage  and  Debye  length  which  is  approximately  valid  for  a 
large  range  of  values  of  these  two  parameters.  The  arguments  leading  to  these 
results  are  given  in  the  following  sections. 


2.  THE  SPHERICAL  CASE 

We  first  write  down  the  two  equations  for  a spherical  collector  in  a conven- 
iently  nondimens ionalized  form  thus:  (See  also  previous  references.) 


J.p2  [l - (1  - 1/P2) .-♦/'» 2-''l 

(1) 

J = C2  «/>3/,2/[o(p)]2  /l  . 

(2) 

1 ^ / 

The  definitions  of  the  variables  appearing  in  Eqs.  (1)  and  (2)  are: 


3.  Suits,  C.  G.,  Ed.  (1961)  The  Collected  Works  of  Irving  Langmuir,  Vol.  3, 

Pergamon  Press. 

4.  Chen,  Francis  F.  (1965)  Plasma  Diagnostic  Techniques.  Huddlestone  and 

Leonard,  Eds,  Academic  Press. 
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•A. 


eV„ 

* = vr2 


c = 


M 2vfir 


A /r 
D'  p 


and 

r = probe  radius 
P 

r = sheath  radius 
s 

Ap  = Debye  length 
Vp  = probe  potential 
e = electronic  charge 

N = ambient  undisturbed  charged  particle  density 
00 

T = (Maxwellian)  temperature  of  ambient  charged  particles 
m - mass  of  collected  particle 
k = Boltzmann  constant 
I = Probe  current 

We  note  that  a is  a known  function  of  p,  so  that  Eqs.  (1)  and  (2)  are  two  equa- 
tions in  the  four  variables  C,  0,  p and  J.  Concerning  Eq.  (1),  the  LMS  equation, 
there  is  in  its  derivation  no  assumption,  aside  from  an  indirect  one  related  to  the 
shape  of  the  potential,  relative  to  the  size  of  the  sheath.  There  is  only  an  assump- 
tion that  a sheath  edge  exists,  beyond  which  the  plasma  properties  assume  their 
undisturbed  values.  Thus,  in  the  absence  of  other  restrictions,  Eq.  (1)  is  expected 
to  be  valid  for  any  sheath  size.  Concerning  Eq.  (2),  the  SCL  equation,  the  limiting 

assumption  is  that  all  particles  which  leave  the  sheath  edge  arrives  at  the  probe. 

2 2 
Now,  when  <t>/(p  - 1)  » 1 in  Eq.  (1),  the  latter  reduces  to  J = p : that  is,  all  the 

particles  leaving  the  sheath  edge  reaches  the  probe.  In  fact,  when  the  more 
2 

severe  condition  0/p  » 1 holds,  the  motion  of  charged  particles  to  the  probe  is 

approximately  radial.  The  key  to  the  simultaneous  validity  of  Eqs.  (1)  and  (2)  is 

thus  the  criterion  0/{p^  - 1)  » 1.  When  0/(p^  - 1)  « 1 on  the  other  hand,  Eq.  (1) 

reduces  to  J = 1+0  which  is  independent  of  p.  We  thus  arrive  at  a situation  in 

which  solving  the  two  equations  simultaneously  for  J as  a function  of  0 and  C yields 

2 2 

results  which  are  valid  both  for  0/(p  - 1)  >>  1 and  for  0/(p  - 1)  « 1.  Since 

o 

0 / (p  - 1)  = 1 is  not  a singular  point  for  either  equation  and  further  since  Eq.  (2) 

2 

is  not  too  bad  an  approximation  at  0/(p  - 1)  = 1,  it  is  not  unreasonable  to  assume 

2 

that  values  of  J obtained  for  0/(p  - 1)  Si  are  also  approximately  correct.  As  far 


9 


1 


as  values  of  p are  concerned,  however,  this  is  clearly  not  the  case;  the  values  of 

2 

p obtained  are  only  accurate  for  </>/(p  - 1)  » 1.  We  recognize  the  heuristic 

* 

nature  of  these  arguments  and  erolore  the  consequences  of  their  assumed  validity. 

To  solve  Eqs.  (1)  and  (2)  we  form: 

Ftp)  = p2  [l  - (1  - 1/p2)  e'^/(P2_l)]  [Q(p)l  2 - C2  *3/2  ^1  + = 0 . (3) 

Using  the  Newton -Raphson  method  Eq.  (3)  is  solved  for  p,  and  J is  then  obtained 

from  either  Eq.  (1)  or  (2).  Details  of  this  solution  are  given  in  the  appendix.  1 

The  results  of  these  calculations  are  shown  in  Figures  1 and  2.  In  Figure  2, 
values  of  p - 1 above  and  to  the  left  of  the  curve  labeled  p =V1  + <t>  are  available, 
but  as  noted  the  calculation  does  not  yield  valid  results  for  p in  this  region.  As 

the  curve  p = >/l  + <j>  is  approached,  the  accuracy  of  the  values  for  p is  reduced,  i 

while  values  away  from  this  curve  will  have  increasingly  better  accuracy.  The 

dashed  line  in  Figure  1 is  for  J = 1 + </>(l  - 1/e)  and  is  hence  the  equivalent  of  the  | 

curve  p = vT  + <t>  in  Figure  2. 

J 

1 


Figure  1.  Nondimensionalized  Current  J vs  Nondimension  - 
alized  Potential  ^ for  Various  Values  of  the  Parameter  C for 
a Spherical  Collector 
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Figure  2.  Nondimensionaiized  Sheath  Size  p vs  Nondimen- 
sionalized  Potential  <j>  for  Various  Values  of  the  Parameter 
C for  a Spherical  Collector 


3.  ASYMPTOTIC  PROPERTIES 


We  investigate  the  dependence  of  J and  p on  <j>  and  C as  <f>  — oo,  and  as  <j>  — 
Note  the  difference  between  this  asymptotic  limit  and  the  limits  commonly  dis- 
cussed in  probe  theory.  * In  the  latter,  <j>  is  held  fixed  while  p — oo  or  0;  here, 
— oo  or  0 while  p varies  as  an  implicit  function  of  The  procedure  followed 
here  is  the  more  useful,  since  the  variables  under  direct  experimental  control 
are  <j>  and  C,  not  p. 

As  <t>  -»  oo,  we  may  assume  that  p — oo  so  that  c?  — (1. II)2/2  p2 ^2  (See 
Appendix,  Eq.  (12)).  Thus,  Eq.  (2)  may  be'  replaced  in  this  limit  by 


J = 


_cV/2 

(i.  nr2  p3/2 


There  are  now  three  possibilities. 

(1)  As  <j>,  p -*  oo,  <|>/(p2  - 1)  -»  oo 

(2)  As  4,  p - oo,  <*/(p2  - 1)  - 0 

(3)  As  0,  p -»  oo,  </>  / (p  ^ - 1)  -»  _k_  * 0,  oo  . 


I 


I 

I 

* 


We  examine  each  of  these  in  turn  to  see  if  th  " can  represent  self-consistent 
solutions  of  Eqs.  (1)  and  (4). 

(1)  In  this  case  4 » p2  - 1,  so  that  J = p2  (Eq.  (1)),  p ' ^2  - C2  4 3 /{l.  1 l)3 

(Eq.  (4)),  4 = 1.  11  p7/3/C4/3,  <t>/(p2  - 1)  as  1.  11  p1/'3/C4/'3  - oo  and  this  is 
self-consistent. 

We  have  for  this  case 

p =C4/V/?/<  l.ll)3/7  (5) 

J = C8/7  *6/7/(l.  ID6/7  . (6) 

(2)  p2  - 1 » Then,  J = 1 + 4 - 4 (Eq.  (D),  p3/2  = C2  41/2/(l.  ll)3^2 

(Eq.  (4)),  4 = (1.  ID3  p3/C4,  4 /(p 2 - 1)  “ (1.  ID3  p/C4  - 00  and  this  is  not 
self  consistent.  2 — 

(3)  4 = k(p2  - 1).  Then,  ,llm  e"^'5  = e'-s  B a constant, 

J = p2(l  - B)  + B (Eq.  (D).  p2(l  - B)  + B = C2/(l.  11)3/2  43/2/p3/2,  and  for  p 

large  this  is,  outside  of  the  factor  (1  - B)  identical  with  case  (1),  which  leads  to 
2 

4/(p  - 1)  — 00  and  is  hence  not  self-consistent. 

Thus,  only  case  (1)  is  self-consistent  and  Eqs.  (5)  and  (6)  must  be  the 
asymptotic  limits  for  4 — « . 

The  two  conditions  necessary  for  Eqs.  (5)  and  (6)  to  be  valid,  namely 
2 

1.  lip  » 1.  64  and  4/(p  - 1)  » 1 are  asymptotically  equivalent  to 

(1.  11)4/7  c4/7  43/7  » l.  64  and  (1.  ll)3^7  4*^7/C3^7  » 1 respectively.  The 
accuracy  of  the  asymptotic  approximation  thus  increases  with  increasing  4,  but 
may  either  increase  or  decrease  with  increasing  C.  This  behavior  is  confirmed 
by  a comparison  of  numerical  and  asymptotic  values.  The  agreement  is  best  at 

C = 1 in  which  case  the  asymptotic  value  is  accurate  to  within  2 percent  at 

4 6 

4 = 10  and  to  within  0.  4 percent  at  4 = 10  for  both  J and  p. 

The  opposite  limit,  4 — 0 can  also  be  investigated  using  Eq.  (2)  since  the 

latter  is  derived  for  charged  particles  leaving  the  emitter  surface  with  nonzero 

initial  velocity,  and  is  hence  valid  for  small  values  of  4-  Although  this  case  is 

not  of  great  interest  practically,  a consistency  analysis  similar  to  that  performed 

for  the  4 — 00  limit  is  easily  carried  out  and  shows  that  the  only  consistent  results 

for  4 — 0 are 


4 /<P  2 - l) 


4*/2 

2s/TT66  C 


and  hence  that 


lim 

4-0 


J = 1 + 4. 


12 


For  fixed  <£,  the  approximation  becomes  better  with  increasing  values  of  C.  Thus, 
for  example,  for  <£  = 0.01,  the  asymptotic  value  of  J is  1.01  and  the  numerical 
values  are  as  follows: 


c 

0.  01 

0.  03 

0.  1 

0.3 

1.0  3 

10  30 

100 

J 

1.  003 

1.  006 

1.  009 

1.  010 

1.010  1.010 

1.010  1.010 

1.  010 



2 

In  this  limit  (</>  - 0)  since  $/(p  - 1)  « 1,  there  is  no  valid  approximation  for 

P- 


4.  THE  CYLINDRICAL  CASE 


The  argument  for  the  cylindrical  case  is  not  as  clearcut  as  is  that  for  the 
sphere.  The  LMS  equation  for  the  cylinder  contains  <j>  in  three  combinations  with 
p and  an  addition  assumption  is  needed  to  complete  the  argument,  namely  > 1. 
Nevertheless,  values  of  <ji  for  which  this  holds  are  the  ones  of  greatest  practical 

interest. 

The  two  applicable  equations  for  the  cylindrical  case  are  (previou-  references) 


J = p erf 


J 


= C 


03/2 

13(p)1  2 


+ 

v/0  / 


where 


J=l[2'LrpeNB 


-1 


(7) 


(8) 


L = probe  length  and  all  other  quantities  are  as  previously  defined  for  the  sphere. 

2 2 2 

Now,  let  0 > 1 , and  to  start  with,  0/(p  - 1)  « 1.  Then,  p 0 /(p  - 1)  a <j>  and 

Eq.  (7)  becomes 


Jap  -2—  \l  — + e^  ( 1 - erf  sfy  ) 

^ ? pi  - 1 


a -H_  + e^  (1  - erf  ) 

Jif 


13 


2 

which  is  independent  of  p.  In  the  opposite  case,  namely  i J>/(p  - 1)  » 1,  we  must 

2 2 

also  have  p <f>/(p  - 1)  » 1 and  hence  JSp,  Thus,  with  the  restriction  <j>  > 1,  the 

same  reasoning  as  applied  to  the  spherical  case  applies  here,  and  the  two  Eqs.(7) 

and  (8)  may  be  solved  simultaneously  to  obtain  J and  p as  functions  of  <j>  and  C. 

2 

Again,  the  solutions  for  p are  only  valid  when  <^/(p  - 1)  » 1. 

To  solve  Eqs.  (7)  and  (8)  we  form 


cy/2 

t3(p)l 1 


(9) 


solve  this  for  p,  and  then  obtain  J from  either  Eq.  (7)  or  (8).  The  details  of  the 
numerical  solution  are  again  given  in  the  appendix. 

The  results  of  these  calculations  are  shown  in  Figures  3 and  4.  It  is  seen 
that  for  the  cylindrical  collector,  the  solutions  for  p have  a domain  of  validity  which 
differs  markedly  from  the  domain  of  validity  of  the  sphere. 

The  asymptotic  properties  of  the  solution  for  the  cylinder  may  be  obtained  by 
considerations  of  self-consistency  in  the  same  manner  as  was  done  for  the  sphere. 

The  results  of  such  considerations  are  again  that  there  is  only  one  self-consistent 

2 

solution  for  the  two  equations  and  for  0 -*  oo,  this  is  that  <f> / (p  - 1)  -*  0,  and 


lim 
<t>  -*  00 


J = 


(10) 


The  asymptotic  expression  for  p may  also  be  formally  written  down  but  is 

2 

meaningless  since  solutions  for  p are  only  valid  when  $/(p  - 1)  » 1,  while  the 

2 

asymptotic  forms  are  only  valid  for  /(p  - 1)  « 1.  The  asymptotic  expression 

for  p is  thus  valid  nowhere. 

Note  the  radical  difference  in  behavior  between  the  spherical  and  cylindrical 
cases:  for  the  sphere. 


lim 
lj>  -»  00 


4/(p2 


1) 


00 


but  for  the  cylinder. 


lim 
d>  -*  oo 


<t>/(p2 


1)  - 0 . 


The  high  voltage  limit  of  the  current  for  the  sphere  is  thus  space  charge  limited, 
while  for  the  cylinder  the  high  voltage  limit  is  orbit  limited. 


Since  for  the  cylinder  we  must  in  any  event  have  <j  ^ 1,  there  can  be  no 
asymptotic  forms  for  <j>  — 0. 

2 

The  conditions  for  which  Eq.  (10)  is  valid,  namely  <j>  » 1,  b/P  « 1 indicate 
that  the  accuracy  of  the  asymptotic  formula  increases  with  increasing  C and  also 
with  increasing  <j>  for  higher  values  of  This  behavior  is  evident  in  Figure  3,  and 
is  also  confirmed  by  the  numbers  below  which  give  percent  deviations  between 
asymptotic  and  numerical  values  of  J for  selected  values  of  C and  <fr. 


c 

3. 

3. 

3. 

3. 

1. 

1. 

1. 

.3 

. 3 

<t> 

102 

103 

104 

105 

103 

104 

105 

104 

105 

percent 
dev  J 

0.  89 

0.  14 

0,00 

0.  00 

4.  23 

1.  34 

0,34 

36.  4 

Appendix  A 

Formulas  for  Numerical  Calculations 

Listed  below  for  reference  purposes  are  the  expressions  for  a,  a',  3,  3'  and 
F'  used  in  the  Newton-Raphson  solutions  of  Eqs.  (3)  and  (9).  (Primes  denote 
derivatives  with  respect  to  p,  and  y = In  p. ) 

(1)  Spherical  Case.  The  values  of  trip)  needed  to  solve  Eq.  (3)  are  tabulated, 
but  we  used  a series'' 

ax  = y + 0.  3 y1 2  + 0.  75  y3  + 0.  143182  74  + 0.  0021609  75  + 0.  00026791  7 6 

(Al> 

2 

for  values  of  p f 10.  0 and  a formula 

o2  = (1.  lip  - 1.64)3/4  (A2) 

For  values  of  p > 10.0.  At  p = 10,  Eq,  (Al)  yields  o2  = 29.06,  Eq.  (A2)  yields 
£*2  = 29.  09,  and  the  tabulated  value  is  a2  - 29. 19.  The  value  of  p is  tested  at  each 
iteration  in  the  Newton-Raphson  process  and  or  a2  chosen  depending  on  whether 
p < 10  or  p > 10.  The  other  expressions  needed  for  the  spherical  case  are: 

1.  Suits,  C.G. , Ed.  (1961)  The  Collected  Works  of  Irving  Langmuir,  Vol.  4, 

Pergamon  Press. 

2.  Suits,  C.G. , Ed.  (1961)  The  Collected  Works  of  Irving  Langmuir,  Vol.  3, 

Pergamon  Press. 


17 


F'  = p2  [l  - (1  - 1/p2)  e'^/(p  _1)]  2 ns' 

+ 2 a2  p [l  - (</,/(p2  - 1)  + 1)  e'^P2'1*] 

«'  = I [l.0  + 0.6r  + 0.  225  y2  + 0.  0572728  y3  + 0.  0108045  y4  + 0.  00160746  y5] 
1 P 

a'2  = 0.8725  (1.  lip  - 1.64)*1/4  . 

2 

(2)  Cylindrical  Case.  For  /3(p)  we  used 

= y + 0.  40  y2  + 0.  09166C7  y3  + 0.  0142424  y 4 + 0.  0016793  y5  + 0.  0001612  y6 

(A3) 

for  p < 16  and 

32  = 4.  6712P(log10  p - 0.  1505)3/2  (A4) 

for  p > 16.  When  p = 16,  Eq.  (A3)  gives  02  = 80.  83,  Eq.  (A4)  gives  02  = 80.  83, 
and  the  tabulated  value  is  81.  203.  The  remaining  expressions  for  the  cylindrical 
case  are: 


F'  = - 2/s/*  _1)  s/^/Vp2- 1 + erf, 


P2-1 


+ C2^3/2  1 + 2J6  2 3, 

\ 4 / r 


0^  = I [l  + 0.  8 y + 0.  275  y2  + 0.  0569697  y3  + 0.  0083965  y4  + 0.  0009672  y5] 


P2  = 2.3356  j-  [0.  65144  (log10  p - 0.  1505)1/2  + (log1Q  p - 0.  1505)3/2] 
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